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Identical boson Hanbury-Brown-Twiss interferometry as applied to relativistic heavy-ion 
collisions is reviewed. Emphasis is placed on the use of hadronic scattering models to 
interpret the physical significance of experimental results. Interferometric studies with 
center-of-mass energies from < 1 GeV/nucleon up to 5500 GeV/nucleon are considered. 



1. Introduction 

1.1. Scope of this Review 

In the present work the apphcation of the Hanbury-Brown-Twiss interferometry 
(HBT) technique to relativistic heavy-ion colhsions (RHI) is reviewed. The em- 
phasis is placed on comparing identical boson HBT experiments ranging in ^/s 
from less than 1 GeV/nucleon to 200 GeV/nucleon with hadronic scattering mod- 
els to attempt to understand what has been learned from these studies over the 
last twenty-five years. Predictions from such a model are also given for future LHC 
Pb+Pb collisions at 5500 GeV/nucleon. Since the literature has been quite rich in 
this field during this period, it has been necessary to be selective both in the ex- 
periments and models which are presented. Thus, while in no way attempting to be 
comprehensive, an attempt has been made to present results in this work which are 
at least representative of the major developments in the field. More comprehensive 
reviews of this field can be found in the literature. ^'^ 

The review is organized as follows: the remainder of Section 1 gives a brief 
discussion of the original work of Hanbury-Brown and Twiss in developing HBT 
followed by the motivation for applying HBT to RHI, Section 2 gives some practical 
information to help understand how HBT is applied to RHI, Section 3 describes 
the hadronic scattering models discussed in this paper and presents comparisons of 
these models with results from two-boson HBT experiments, and Section 4 gives a 
summary. 
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1.2. Origins of HBT: Hanbury-Brown and Twiss experiments 

About fifty years ago, Hanbury-Brown and Twiss first suggested, and then proved in 
a table-top experiment, that photon pairs exhibit a second-order interference effect 
if detected simultaneously in two detectors.^ They applied this technique, which 
we now call HBT when applied with any type of boson,** to the measurement of 
the angular diameter of stars using pairs of photomultiplier tubes to detect optical 
photons^ and pairs of radiotelescope dishes to detect longer wavelength photons.^ 
Their measurement of the optical angular diameter of the star Sirius located in 
the constellation Canis Major, serves as a good example of the HBT method.^ A 
schematic layout which helps to demonstrate the principle of their method is shown 
in Figure 1 . Sirius is shown emitting two photons with wavevectors and kj from 
points Xj and Xj, respectively, which are detected in two photomultiplier tubes 
located at positions ai and a2 on the Earth. 



PM 
tubes 




Fig. 1. Schematic diagram of the method used to measure the angular diameter of Sirius by 
Hanbury-Brown and Twiss in 1956. 

Assiiming that photons arc emitted incoherently from the star at each position 
and the photons are emitted as plane waves and taking a time-independent picture, 
the wavefunction to detect in coincidence the two photons in the two detectors on 
Earth, ^'(ki, k^; Xj, Xj), is 

*(kj, kj; Xj, Xj) = 6[exp(ikj • (ai - Xj)) exp(ikj • (a2 - Xj)) 

-I- exp(ikj • (a2 - Xj)) exp(ikj • (ai - Xj))] (1) 

where the second term arises due to the path ambiguity of detecting bosons (as 

shown in Figure 1), and where 6 is a normalization constant. The probability to 
detect the two photons, Pij, is just the square of ^P, i.e. and thus given by 

Pij ( Ak, d) = 6^ [1 -h cos(Ak • d)] (2) 

where Ak = k^ — kj and d = ai — a2. The separation between the two detectors, 
d, is called the baseline. It is seen in Eq. (2) that a) the cos term in P^j is a result 
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of the path ambiguity term in Eq. (1), and b) Pij does not depend on the photon 
emission positions at the star but only on the differences between the wavevectors 
and the detector positions. One can write Eq. (2) in a more convenient form with 
the approximations ki ~ kj = k, Ak • d ~| Ak |, | Ak |~ Oijk = 2Tr9ij/X, taking 
b=l, and defining the correlator, Cij{d) = Pij — 1, resulting in 

Cy(rf) = cos(27r%d/A) (3) 

where, as seen in Figure 1, Oij is the angular diameter seen on Earth between the 
points on the star i and j and A is the wavelength of the photons. Cij {d) is propor- 
tional to the coincidence signal produced in the photomultipliers for the photons 
coming from these two points on the star, but in practice, the signal measured in 
the electronics is the sum over all of the photon emission point pairs from the star. 
If N is defined as the number of photon source points making up a star of brightness 
(intensity) / and radius r, to get the detected signal from the entire star, C{d), Cij 
is summed over all of the unique pairs of photon emission points 

N N 

C{d) = e ^ (d) = e ^ cos(27r%(i/A) (4) 

i>j i>j 

where e is a conversion factor to get the detector signal. If one looks at cases where 

the arguments of the cos terms in Eq. (4) are small such that the cos terms arc not 
far from unity, e.g. for small baselines, Eq. (4) can approximately be expressed as 

C{d) ~ e— cos(27r6lirf/A) (5) 

where is the average of the % . 

One is now in the position to measure the average angular size of a star by 

measuring C{d) as a function of d and fitting these data with Eq. (5) to extract 
9i. Figure 2 shows measurements of C{d) for Sirius by Hanbury- Brown and Twiss. 
Also plotted are fits to the data using Eq. (5) (taking A = 550 nm) and a full 
calculation including integrations over the photon wavelength spectrum and photon 
source distribution and detection efficiency effects by Hanbury-Brown and Twiss. ^ 
As seen, 9i extracted from the fit to the data by Eq. (5) agrees with the angular 
diameter extracted from the fit by the full calculation within a factor of two. 

Another example of applying the HBT technique to measure angular size is for 
the binary star Alpha Centauri, located in the constellation Centauris seen in the 
southern hemisphere. A schematic diagram of the geometry of the measurement is 
shown in Figure 3. For the orientation of the two stars shown, there are two angular 
size scales that enter this measurement: the angular sizes of the individual stars 
(assumed to be identical), 9i, and the angular size between the two stars, 02- What 
angular size will the HBT technique measure? To answer this question, calculate 
Cbinary{d) using Eq. (4) assuming that the binary star has the total number of 
photon source points N distributed equally between the two stars and setting e = 1 
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Fig. 2. Measurement of C{d)/C{0) vs. d for Sirius by Hanbury-Brown and Twiss. Also shown are 
fits to the measurements using Equation 4 and a full calculation. 

for simplicity, giving 

N/2 N/2 
Cbinaryjd) = 2 COs(27r%(j/A) + ^ COs{2TT9ijd/ X) (6) 
i>j i{A),j(B) 

where the first sum accoimts for the two identical sums over the individual stars 
in the binary, and the second sum accounts for the sum over photon source pairs 
between the two stars, i.e. the index i associated with star A and j with the star 
B. One can now use the same approximation as before in going from Eq. (4) to 
Eq. (5), i.e. the argument of the cosine being small, to express the first sum in Eq. 
(6) in terms of 9i. This will not be the case for the second sum, since for Alpha 
Centauri the separation between the stars, L, is much larger than the radius of 
the individual stars, r (_L ~ 11 A.U. and r ^ solar radius), resulting in 02 » Oi. 
Thus, the cos terms in the second sum will oscillate wildly, resulting in that sum 
vanishing, and Eq. (6) becomes 

Cbinary {d) = —COs{2neid/X). (7) 

Comparing Eq. (7) with Eq. (5), it is seen that the binary star gives half the detected 
HBT signal expected from a star with N photon source points, and its angular size 
is characterized by the size of the individual stars making up the binary. When 
Hanbury-Brown and Twiss performed such a measurement on Alpha Centauri, this 
is indeed what they observed.^ 



I 



I 
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Fig. 3. Schematic diagram of the geometry for the measurement of the angular size of the binary 
star, Alpha Centauri. 

This example of the use of this technique on a binary star shows that the par- 
ticular boson source distribution of the system measured has an important impact 

on how results from HBT arc interpreted. It will be shown later that similar effects 
are seen when applying this technique to two-pion measurements of RHI. 

2. Applying HBT to Relativistic Heavy-Ion Collisions 

Relativistic heavy-ion collisions are similar to stars in that they emit bosons from 
a finite-sized region (i.e. the interaction region). They differ from stars in that 1) 
the bosons are predominantly pions, 2) other bosons such as kaons are emitted, 
3) the size scale is much smaller (~ 10~^^ m), and 4) the lifetime of the boson- 
emitting source is short 10^^'' s), introducing time as an important variable. 
Nevertheless, one could imagine applying HBT to RHI to directly measure the size, 
and perhaps the lifetime, of these collisions and thus get information on their space- 
time dynamics. The first measurement of the pion-emitting source using two-pion 
interferometry was carried out by G. Goldhaber, S. Goldhabcr, W. Lee. and A. 
Pais in 1960 with ^ = 2.1 GeV proton-|-antiproton collisions from the Bevatron 
(they referred to this method as Bose-Einstein correlations and others have called 
this the GGLP effect).^ In this section a qualitative derivation is given for two-pion 
Bose-Einstein/HBT interferometry and a few practical considerations for carrying 
out experimental RHI HBT studies are given. 

2.1. Qualitative derivation of two-pion HBT in RHI 

A qualitative derivation of two-pion HBT interferometry is given below. More for- 
mal derivations can be found in the literature. ^'^ As will be seen, its derivation 



starB 
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will closely parallel that given earlier for two-photon HBT. A schematic diagram of 
the geometry of a two-pion interferometric measurement is given in Figure 4. The 

interaction region in a RHI collision is shown emitting two pions from points and 
Xj , which are detected with wavevectors kj and kj in two "pion detectors" located 
at positions ai and a2, respectively, in the experimental hall. 




Fig. 4. Schematic diagram of a two-pion interferometric measurement in which the pions are 
produced in a RHI collision. 

Assuming that pions are emitted incoherently from the interaction region at each 
position and the pions are emitted as plane waves and taking a time-independent 
picture, the wavefunction to detect in coincidence the two pions with these wavevec- 
tors in the two detectors, 4'(kj,kj;Xj,Xj), is 

*(kj, k^; Xj, Xj) = 6[exp(ikj • (ai - Xj)) exp{i\s.j ■ (a2 - x^)) 

-I- exp(ikj • (ai - x^)) exp(ikj • (a2 - Xj))] (8) 

where the second term arises due to the path ambiguity of detecting bosons (as 
shown in Figure 4), and where 5 is a normalization constant. This expression is 
similar to Eq. (1) above used for photon pairs from a star, but not identical. The 
difference is that indices are switched in the second term since we assume the pion 
wavevector (momentum) is now associated with a particular detector, rather than 
with an emission point as in the case of a photon from the star. The ambiguity in 
this case is in which emission point the pion originates from, not which detector 
records a photon. The probability to detect the two pions, Pij, is just the square of 
^, i.e. ^*^, and thus given by 

Pij ( Ak, r) = 6^ [1 + cos(Ak • r)] (9) 

where Ak = k^ — kj and r = x^ — Xj. The separation between the two emission 
points, r, is related to the size of the interaction region or "pion source." It is seen 
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in Eq. (9) that a) the cos term in Pij is a result of the path ambiguity term in Eq. 
(8), and b) Ptj does not depend on the detector positions but only on the differences 

between the pion momenta and the emission positions. Note that the time variable 
can also be introduced into this formalism by considering k and x as four vectors. 

As in the discussion of HBT with stars in the last section, the interaction region 
can emit pions from an extended region of points in space which the detectors smn 
over. Defining the pion source distribution function as a contimiOTis function of the 
position of a source point x, p(x), one can integrate the interaction region over all 
pairs of source points using Eq. (9), 



where C(Ak) is called the correlation function. Assuming for simplicity that the 
pion source is distributed as a gaussian of width R, i.e. p{r) ~ exp(— r^/2i?^), and 
substituting into Eq.(lO) the result is, 



Figure 5. shows a two-pion HBT measurement in central RHIC Au+Au collisions 
from the STAR experiment, which detects charged particles in a large-acceptance 
Time Projection Chamber immersed in a magnetic field. The one-dimensional 
two-pion correlation function (extracted from the two-pion coincidence count rate, 
see later) evaluated in the invariant frame of the pion source is plotted versus Ak. 
Only the data points lower than HAk < 50 MeV/c are plotted. A version of Eq. 
(9) assuming Ak ± r including an arbitrary parameter A is fitted to the data to 
extract r and A (dashed line). A version of Eq. (11) with the same parameters is 
also fitted (sohd line). As seen, the gaussian function fits the data better, but both 
functions give about the same fitted values of r ~ 8 fm and A ~ 0.5 — 0.6. 

What is the interpretation of these results for r and A? Naively, one would 
expect the size of the Au nucleus, 1.2A^^^ ~ 7 fm, to set the scale for the size of 
the pion emitting source in these collisions, and this is seen to be the case for r. 
For A one might expect this to be unity for the idealized derivation of Eq. (9), yet 
the value extracted in this case is significantly smaller. This is suggestive of the 
reduced HBT signal seen in binary stars discussed above, for which boson sources 
of distinctly different sizes are intermixed. 

How to interpret the parameters from such fits to data in RHl HBT analyses 
will be an important theme for the rest of this paper. As will be shown, models of 
the interaction will be used as tools to interpret experimental HBT results. 

2.2. Measuring the Correlation Function in experiments 

As shown in Figure 4., in carrying out an experimental HBT study one measures the 

two-boson coincident count rate along with the one-boson count rates for reference. 
The experimental two-boson correlation function for bosons binned in momenta ki 




(10) 



C(Ak) ~ l+exp[-{AkfR^]. 



(11) 



February 5, 2008 15:29 WSPC/INSTRUCTION FILE hbtreviewf 



8 Thomas J. Humanic 




10 20 30 40 



Ak*(h/27t) (MeV/c) 



Fig. 5. Two-pion HBT in RHIC Au+Au collisions from the STAR experiment with fits. Only 
points for HAk < 50 MeV/c are plotted. 



and k2, C(ki,k2), is constructed from the coincident countrate, 7V2(ki,k2) and 
one-boson countrate, Afi(k), as 

C(k.,M.a,k.,k.,^^^^|^, (12) 

where a(ki,k2) is a correction factor for non-HBT effects which may be present 
in the experiment. Typically,the largest contributions to Q:(ki,k2) occur in the 
correction for the boson detection efficiency and in correcting for final-state boson- 
boson Coulomb repulsion.^ It is usually convenient to express the six- dimensional 
C(ki, k2) in terms of the four- vector momentum difference, Q = \ki — k2\ hy sum- 
ming Eq.(12) over momentum difference, 

where A{Q) represents the "real" coincident two-boson countrate, B{Q) the "back- 
ground" two-boson countrate composed of products of the one-boson countrates 
and e{Q) the correction factor, all expressed in Q. ^ In practice, B{Q) is the mixed 
event distribution, which is computed by taking single bosons from separate events.^ 



I 



I 
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2.3. Parameterizing the Correlation Function: extracting boson 
source parameters 

The experimental two-boson correlation function is formed using Eq.(13) and a 
gaussian model for the boson source distribution is normally fitted to it to extract 
radius parameters. A standard parameterization of C(Q) is obtained by assuming a 
gaussian space-time distribution of freeze-out points, p{r, t), in terms of the variables 
fTsidei which points in the direction of the sum of the two boson momenta in the 
transverse plane, rxout, which points perpendicular to rrside in the transverse plane, 
the longitudinal variable along the beam direction, z, and time, t: 

p(r,i) = Fexp(--^;2^--^J^-— ^-^) (14) 

^^Tside "^^Tout "^^Long 

where I^Tside ^ transverse sideward radius parameter, R'xout ^ transverse out- 
ward radius parameter, R'l^ong ^ longitudinal radius parameter, r is a lifetime 
parameter, and F is a normalization constant. From this distribution function 
the following parameterization of C(Q) can be obtained 

C{QTside, Qrout, Qhong) = 
G[l -h y^exJ£>{ — QTside^side ~ Qrout^Tout ~ QLong^Long)] (l^) 

where Q has been broken up into the two transverse and one longitudinal com- 
ponents, G is a normalization constant, and A is the usual empirical parameter 
added to help in the fitting of Eq. (15) to the actual correlation function (A = 1 in 
the ideal case). The radius parameters in Eq. (15) are related to those in Eq.(14) 
in, for example, the LCMS frame (longitudinally comoving system in which the 
longitudinal boson pair momentum vanishes) as follows; 

^Tside ~ ^Tside 

Rrout = R'tou' + (16) 

^Long ^Long 

where Pt is the transverse velocity of the boson pair. In this parameterization as 
seen from Eq. (16) Rxout contains information about both the transverse size and 
the lifetime of the source. Note that Eq.(15) follows from Eq.(14) under the as- 
sumption of a geometrically static boson source. In a realistic heavy-ion collision 
the source will not be static and may have position-momentum correlations and 
other effects that could make the source parameters defined above depend on the 
boson pair momentum . In the present application, Eq. (15) is fitted to the exper- 
imental correlation function to extract the radius parameters R.Tside, Rtoui, and 
Rhong and A. Figure 6 shows experimental 7r~ — 7r~ correlation functions with and 
without final-state Coulomb corrections included from the Relativistic Heavy Ion 
Collider (RHIC) STAR experiments^ for ^ = 130 GeV/nucleon Au-|-Au collisions 
along with gaussian fits. 
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Fig. 6. TT —IT HBT in RHIC Au+Au collisions from the STAR experiment with gaussian fits. 
The Coulomb-corrected data are shown as filled circles and the uncorrected data as open circles 
(from Ref.[10]). 



3. Comparing hadronic scattering models with two-boson HBT 
experiments from RHI collisions 

Relativistic heavy-ion collisions provide a means of creating matter in a hot and 
dense state which might shed light on the behavior of matter under these extreme 
conditions, such as the possibility of producing a phase transition to quark matter. 
It is generally agreed that the most extreme conditions exist in the initial state of 
the heavy-ion collision, roughly defined as occurring just after the projectile and 
target nuclei pass through each other. Eventually the interaction region hadronizes 
into a large number of mesons and baryons (with non-hadronic particles such as 
photons, electrons, and muons also being produced) and then expands to its final 
state. During the expansion stage the hadrons strongly scatter with each other, this 
process being called rescattering. The final state of the collision can be thought of 
as the state for which rescattering ceases among all remaining (final) hadrons. It is 
often convenient to define the "freezeout" point of a final hadron as the position, 
time, energy, and momentum the particle had when it stopped rescattering. Thus, 
one can define more precisely the final state of the collision as the collection of 
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the freezeout points of all of the final hadrons in an 8- dimensional phase space 
(4 dimensions for space-time and 4 dimensions for momentum-energy). The term 
freezeout will be used to represent the final state of the collision defined in this 
way. In principle, the properties of the collision at freezeout are directly accessible 
by measurement. In practice, one directly measures the freezeout momenta and 
energies of the particles using, for example, magnetic spectrometers and indirectly 
measures the freezeout space-time from the freezeout momenta and energies using 
the method of two-boson HBT interferometry. 

However, the main motivation to study relativistic heavy-ion collisions is to 
obtain information about the initial, extreme state of the collision. Some directly 
measurable non- hadronic probes such as direct photons, and electron and muon 
pairs have been predicted to be sensitive to certain features of the initial state. 
It is equally important to study the final hadrons from the collision since they 
should contain information about the bulk properties of the initial state, i.e. the 
temperature and energy density achieved in the collision. The difficulty in using 
the hadrons to extract this information is that the rescattering process masks the 
initial space-time and momentum-energy information by random scattering and 
thus there is no simple connection between the freezeout information obtained in 
experiments and the initial state. 

In this paper, the method employed to approach this problem is to use hadronic 
scattering calculations to disentangle the rescattering effects from the hadronization 
process. The strategy will be to take a simple model for hadronization and propagate 
the initial hadrons via rescattering to freezeout, adjusting the parameters of this 
model to sec if a parameter set can be found where; the freezeout obscrvablcs from 
the calculation agree with those measured in experiments. Within the context of 
such a model, this parameter set thus describes the state of the collision before 
rescattering, putting one a step closer in time to the initial state. The advantage 
of using a simple hadronization model is that the number of parameters to be 
adjusted is minimized, increasing the chances that the extracted parameter set is 
unique. The disadvantage is that some physics of the hadronization will be left out, 
so the physical interpretation of these parameters may be complicated. Descriptions 
of two such hadronic scattering models are presented below. Results of these models 
are then compared with RHI HBT experiments. 

3.1. Intranuclear Cascade Model for y/s < 1 GeV/nucleon: INCM 

Intranuclear cascade models (INCM) have been used extensively to understand gen- 
eral features of RHI collisions for bombarding energies of ^/s < 1 GeV/nucleon.^^ 
The basic assumption in these models is that a RHI collision can be viewed as a su- 
perposition of nucleon-nucleon interactions whose trajectories between interactions 
are described classically while the interactions are determined by experimental scat- 
tering cross sections. The CASCADE code by Cugnon-'^^ is a version of a INCM. In 
this version, A isobars are included which serve as the mechanism through which 
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pions are produced and rescattered throughout the duration of the coUision. The 
CASCADE code is isospin averaged such that there is only one type of pion, nu- 
cleon, and A in the calculation. A pion is defined as freezing out from the calculation 
when it no longer scatters and the time and position of its parent A, along with 
the vector momentum of the pion, are recorded. 

Two-pion HBT predictions are made with this recorded information by weight- 
ing pairs of pions with bose symmetrization using a four-vector version of Eq. (9) 
(so that time effects are included) and then binning the weighted pion pairs simi- 
lar to what is done in experiments using Eq. (13)to form the two-pion correlation 
function. Predictions for pion source parameters are then obtained by fitting a 
gaussian model to this Monte-Carlo generated correlation function. The particular 
pion source model used in this study is similar to that represented by Eq. (14) but 
with a spherical gaussian source, i.e. p{r,t) exp(— "72), giving 

C(Q, Qo) = 1 + Aexp(-^ - ^) (17) 

where the R and r parameters differ by a factor of -\/2 from the similar quantities 
defined in Eq. (14), the A parameter has the same meaning as before, Q = |ki — k2|, 
and Qo = | i?i — i?2 | , where E is the total relativistic energy of a pion. Figure 
7 shows a two-pion correlation function projected onto the Q variable generated 
from this procedure using CASCADE for the reaction 1.5 CcV/nuclcon (lab frame) 
^^4o_j_^j,4o -^yj^jj impact parameter 6 = 2 fm and minimum pion momentum kmm — 
50 MeV/c. The projected fit of Eq. (17) is also shown. The source parameters 
extracted from the model in this case are R = 3.6 ± 0.1 fm, r = 3.2 ± 0.5 fm/c, 
and A = 0.94 ± 0.05. These parameters are shown in the larger context of their 
dependence on kmin and b in Figure 8. As seen, for increasing kmin with fixed 5, 
R and r are seen to decrease somewhat whereas A stays constant at around unity. 
For increasing h and fixed kmin, all three parameters are seen to decrease. 

The mechanisms for these dependencies in the model are seen in Figures 9 and 
10. Figure 9 shows that the freezeout time distribution for larger kmin is shifted 
backward in time compared with smaller kmin ■ Since the pion source from the model 
is found to be smaller at earlier times, this explains the smaller R and t for larger 
kmin'- higher pion momentum probes earlier stages of the collision. Figure 10 shows 
the effect of impact parameter on the size and shape of the pion source. For small 
b the source projected on the x — z plane is spherically symmetric and large, while 
for larger b the source is smaller and extended along the z axis. Thus one expects 
the decrease in R and r for larger b, and sees that the decrease of A is due to the 
source no longer looking like a perfect gaussian. 

3.2. Comparing the IN CM with Bevalac experiments 

In the last section HBT predictions based on the CASCADE INCM were presented 
and analyzed. How seriously to take these results depends on whether or not this 
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1.5 A-GeV'^^Ar + '^^Ar , 
b = 2 FM 




q (MeV/C) 

Fig. 7. Projected two-pion correlation function from CASCADE for the reaction 1.5 GeV/nucleon 
Ar*0 + Ar^O, 6 = 2 fm, kmi„ = 50 MeV/c (From Ref.[17]). 



model gives predictions which agree with experiment. Figure 11 shows a comparison 
of pion source parameters extracted from the model with those from several LBL Be- 
valac HBT experiments: (a) 1.8 GeV/nucleon Ar'^°+KCl and Ne'^°+NaF,^^ (b) 1.5 
GeV/nuclcon Ar-^^ + KCl,^^'"^" and (c) 1.2 GeV/nucleon Ar'^°+KCl.^° The experi- 
ments used two different types of detectors to measure the pion momenta, a narrow 
acceptance magnetic spectrometer^^ and a large acceptance streamer chamber. ^^'^^ 
For the spectrometer, both tt+tt^ and tt^tt^ pairs were used. In order to simulate 
the pion acceptances used in the two experiments, pions used in the model predic- 
tions were selected to be only those which fell in the experimental acceptances. 

As seen in Figure 11, the model predictions, labeled as "Cascade" , mostly agree 
with the experimental source parameters within la of the experimental error bars. 
The largest disagreement is seen for the spectrometer A parameters, the predictions 
being close to unity and the measurements being in the 0.5 — 0.7 range, whereas 
from the streamer chamber experiment A is measured to be close to unity as with 
the predictions. 

From the overall good agreement between model and experiments, one can con- 
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Fig. 8. Dependence of pion source parameters on fcmin and b for the reaction 1.5 GeV/nucleon 
_|_ y5^j.40 "pj^g jjjjgg are drawn to guide the eye (From Ref.[17]). 



elude that the method of symmetrizing the pions produced in the intranuclear 
cascade hadronic scattering model to make HBT predictions can be a valuable tool 
in understanding what HBT is measuring at the "microscopic level" in these lower 
energy RHl collisions. This motivates one to apply this same method to hadronic 
scattering models valid for higher energy RHI collisions. This is done in the next 
section. 



3.3. Hadronic Rescattering Model for y/s > 1 GeV/nucleon: HRM 

For higher energy RHI collisions, i.e. ^/s > 1 GcV/nuclcon, a different hadronic scat- 
tering approach is taken to simulate the collision as the basis for HBT predictions: 
a rescattering calculation is used to disentangle the rescattering effects from the 
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1.5 A-GeV ''^Ar + ''^Ar , 
b = 2 FM 




Time (FM/C) 

Fig. 9. Comparison between final pion freezeout time distributions for kmin 

= 50 and 300 MeV/c, 

for the reaction 1.5 GeV/nucleon Ar"'" + Ar'"' (From Ref.[17]). 

hadronization process. As mentioned earlier, the strategy will be to take a simple 
model for hadronization and propagate the initial hadrons via rescattering to freeze 
out, adjusting the parameters of this model to see if a parameter set can be found 
where the freeze-out observables from the calculation agree with those measured 
in experiments. Descriptions of both the hadronization model and the rescattering 
calculation used are presented below. Results of applying this approach for the BNL 
Alternating Gradient Synchrotron (AGS), CERN Super Proton Synchrotron (SPS), 
BNL Relativistic Heavy fon Collider (RHIC), and CERN Large Hadron Collider 
(LHC) energy collisions are then shown. 

3.3.1. HRM description 

A brief description of the hadron rescattering model (HRM) calculational method 
is given below. A more detailed description is given elsewhere. ^^'^^ Rescattering 
is simulated with a semi-classical Monte Carlo calculation which assumes strong 
binary collisions between hadrons. The Monte Carlo calculation is carried out in 
three stages: 1) initialization and hadronization, 2) rescattering and freeze out, and 
3) calculation of experimental observables. Relativistic kinematics is used through- 
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Fig. 10. Timc-intcgrated spacial pion source distribution projected onto the reaction plane for 
6 = and 6 fm, for the reaction 1.5 GeV/nucleon Ar^*-^ + Ar*" (From R,ef.[17]). 



out. All calculations are made to simulate either AGS, SPS, RHIC, or LHC energy 
collisions in order to compare with present (or future) experimental results 

The hadronization model employs simple parameterizations to describe the ini- 
tial momenta and space-time of the hadrons similar to that used by Herrmann 
and Bertsch.^^ The initial momenta are assumed to follow a thermal transverse 
(perpendicular to the beam direction) momentum distribution for all particles, 

{l/mT)dN/dmT = CmT/[exp (mr/T) ± 1] (18) 

where rriT = \/ Vt^ + "^o^ is the transverse mass, pT is the transverse momentum, 
mo is the particle rest mass, C is a normalization constant, and T is the initial 
"temperature" of the system, and a gaussian rapidity distribution for mesons, 

dNjdy = De^p [~{y - yof/{2ay^)] (19) 

where y = 0.5\n[{E + pz)/{E — p^)] is the rapidity, E is the particle energy, 
is the longitudinal (along the beam direction) momentum, D is a normalization 
constant, yo is the central rapidity value (mid-rapidity), and ay is the rapidity width. 
Two rapidity distributions for baryons have been tried: 1) flat and then falling off 
near beam rapidity and 2) peaked at central rapidity and falling off until beam 
rapidity. Both baryon distributions give about the same results. The initial space- 
time of the hadrons for 6 = fm (i.e. zero impact parameter or central collisions) is 
parameterized as having cylindrical symmetry with respect to the beam axis. The 
transverse particle density dependence is assumed to be that of a projected uniform 
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Fig. 11. Comparison between pion source parameters from CASCADE predictions and pion inter- 
ferometry measurements (a) 1.8 GeV/nucleon Ar"^^ + KCl and Ne^^ -\-NaF, (b) 1.5 GoV/nucloon 
Ar^O + KCl, and (c) 1.2 GeV/nucleon Ar^O + KCl (From Ref.[17]). 



sphere of radius equal to the projectile radius, R {R ~ tqA^^^, where tq = 1.12 
fm and A is the atomic mass number of the projectile). For 6 > (non-central 
collisions) the transverse particle density is that of overlapping projected spheres 
whose centers are separated by a distance b. The longitudinal particle hadronization 
position (zhad) and time {thad) are determined by the relativistic equations,^"* 

Zhad^ Thadsvahy (20) 

thad = ThadCOshy 

where y is the particle rapidity and r^ad is the hadronization proper time. Thus, 
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Table 1. Hadronization model parameters used in the HRM for various systems. 



System 


Thad (fm/c) 


Tinit (MeV) 




SPS Pb+Pb 


1 


270 


1.2 


RHIC Au+Au 


1 


300 


2.4 


LHC Pb+Pb 


1 


500 


4.2 



apart from particle multiplicities, the hadronization model has three free param- 
eters to extract from experiment: <Ty, T and Thad- The hadrons included in the 
calculation are pions, kaons, nucleons and lambdas (tt, K, N, and A), and the p, uj, 
T], T]', (f), A, and K* resonances. For simplicity, the calculation is isospin averaged 
(e.g. no distinction is made among a 7r+, 7r°, and tt"). Resonances are present at 
hadronization and also can be produced as a result of rescattering. Initial resonance 
multiplicity fractions are taken from Herrmann and Bertsch,^^ who extracted re- 
sults from the HELIOS experiment. The initial resonance fractions used in the 
present calculations are: tj/tt = 0.05, p/n = 0.1, p/iv = 3, (f)/{p + uo) = 0.12, 
i /r] = K*/uj = 1 and, for simplicity, A/iV = 0. Note that the AGS Si + Au and 
SPS S + S and S + Pb calculations were done with an earlier version of HRM and 
differ from the description above in that 1) an initial cylinder was used instead of 
Eq. (20) of length z = ±1 fm and 2) no resonances were included. ^^'^^ 

The second stage in the calculation is rescattering which finishes with the freeze 
out and decay of all particles. Starting from the initial stage {t = fm/c), the 
positions of all particles arc allowed to evolve in time in small time steps (dt = 0.1 
fm/c) according to their initial momenta. At each time step each particle is checked 
to see a) if it decays, and b) if it is sufficiently close to another particle to scatter 
with it. Isospin-averaged s-wave and p-wave cross sections for meson scattering are 
obtained from Prakash et al..^* The calculation is carried out to 100 fm/c, although 
most of the rescattering finishes by about 30 fm/c. 

Calculations are carried out assuming initial parameter values and particle mul- 
tiplicities for each type of particle. In the last stage of the calculation, the freeze-out 
and decay momenta and space-times are used to produce observables such as pion, 
kaon, and nucleon multiplicities and transverse momentum and rapidity distribu- 
tions. The values of the initial parameters of the calculation and multiplicities are 
constrained to give observables which agree with available measured hadronic ob- 
servables. As a cross-check on this, the total kinetic energy from the calculation is 
determined and compared with the collision center of mass energy to see that it is in 
reasonable agreement. Particle multiplicities were estimated from charged hadron 
mtiltiplicity measurements and models. The hadronization model parameters used 
for various systems are shown in Table 1. It is interesting to note that it is desirable 
the same value of Thad for all three very disparate in energy systems. 

Figures 12-18 show some results from HRM to give a feeling for the information 
one can obtain from this model. 
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Figures 12 and 13 show the effects of momentum cuts, both in magnitude and 
direction, on the transverse and longitudinal pion source dimensions, respectively. 

Figure 12 shows distributions of pion freeze out positions for central collision SPS 
S-\-Pb collisions projected onto the transverse plane. In (a) and (c) a low momentum 
cut on the pions of < < 100 MeV/c is made, while for (b) and (d) a high pion 
momentum cut of 500 < pr < 600 MeV/c is made. Also, in (c) and (d) only pions 
with azimuthal direction in the range 165° < (j) < 195° (as indicated by the arrows) 
are shown. Comparing (a) and (b), one sees that higher momentum pions tend to 
be concentrated at larger radius compared with the lower momentum pions which 
tend to be imiformly distributed in a spherical volume. As seen in (c) and (d), if 
on top of these momentum cuts on magnitude cuts on direction are also made, the 
geometry of the low momentum pion source is not significantly changed (aside from 
fewer pions satisfying both cuts), whereas the higher momentum pions are seen to 
be more directional and only pions in the vicinity of the angiilar cut are present, 
making the pions source look smaller in size. Figure 13 shows a similar effect along 
the "light cone" of cutting on the direction of pions normal to the longitudinal 
direction, seen in the lower plot of the figure. Since HBT tends to pick out pion 
pairs with small momentum difference, the distributions in Figure 12(c) and (d) 
and Figure 13 give a qualitative indication of how HBT will view these cases. 

One can see the dramatic effect of rcscattcring on the freezeout times of particles 
in the HRM calculation in Figure 14, which shows the freezeout (and decay) time 
distributions for pions, kaons, and nucleons for SPS Pb+Pb collisions (a) from a 
full rescattering calculation, and (b) with rescattering turned off. With rescattering 
the distributions for different particles arc peaked at different freezeout times, the 
kaons peaking earliest and the nucleons the latest. Without rescattering the only 
feature to be seen is the exponential decay of the initial resonances producing pions 
and kaons. 

Figures 15 and 16 show the effects of rescattering from the HRM model on 
the transverse mass distributions for SPS Pb+Pb collisions for pions, kaons, and 
nucleons. Figure 15 shows these distributions at freeze out along with exponential 
fits from which slope parameters are extracted using 

= aexp( — -f) (21) 



mT drtiT T 



where my = ^/pr^ + rn^ is the transverse mass, m is the particle rest mass, T 
is the slope parameter, and a is a normalization constant. Figure 16 shows how 
these slope parameters evolve in time during the calculation. At f = fm/c all of 
the particle types begin with a common temperature (as seen in Table 1) and as 
time proceeds the rcscattcring drives the separation of the slope parameters until 
freezeout at which point they evolve into the experimental values from the NA44 
experiment'' (as discussed above, Tinu is adjusted to give the overall scale of the 
experimental slope parameters). Also shown in Figure 16 is the time evolution for 
the pion slope parameter for the "pion gas" case, i.e. without nucleons and kaons 
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Fig. 12. Distribution of pion freeze out positions for SPS S + Ph collisions projected onto the 
transverse plane, (a) and (c): < < 100 MeV/c pions, (b) and (d): 500 < pr < 600 MeV/c 
pions. In (c) and (d) only pions with azimuthal direction in the range 165° < (j> < 195° (as 
indicated by the arrows) are shown (From Ref. [26]). 

included. As seen, the initial temperature required in this case is much lower than 
for the full calculation, showing the significant effect rescattering with the nucleons, 
and to a lesser extent kaons, has on the pions. 

Figure 17 shows projections of the three-dimensional two-pion correlation func- 
tion for SPS Pb+Pb onto the Qrside, Qtoui, and Qlouq axes from HRM. A fit to Eq. 
(15) is also shown. As seen, applying the boson symmetrization method described 
above with the HRM results in correlation functions which are well described by the 
gaussian model represented in Eq. (15). The two extracted transverse radius param- 
eters, Rxside and Rxout, are seen to be comparable in size around 5 — 6 fm, while 
Rlouq is significantly larger, reflecting the different dynamics of the pion source 
along the beam direction compared with transverse to the beam direction. Looking 
at the resulting A parameter extracted in this case of 0.634 ± 0.019, it is seen to be 
significantly smaller than unity which was obtained in the INCM discussed above. 
The explanation for this in the HRM model is the presence of long-lived resonances 
such as rj and 77'. Pions produced from these resonances come from a much larger 
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Z (fm) 

Fig. 13. Distributions of pion freeze out positions for SPS S + Pb collisions for time versus z. The 
upper plot has no cuts applied. The lower plot only shows positions of pions with radial direction 
in the range 75° < 8 < 105° in the cm. frame (corresponding to a cm. rapidity window of ±0.25) 
(From Ref.[26l). 

source that the directly produced pions, resuhing in an overall suppression of the 
correlation function which is manifest in a smaller overall A parameter for the colli- 
sion. It is interesting to compare this effect with the effect seen by Hanbury-Brown 
and Twiss when observing Alpha Centauri, i.e, Eq. (7). The two cases turn out to be 
analogous since for Alpha Centauri one has two photon sources, i.e. stars, separated 
by a large distance which dilutes the HBT effect, whereas in RHI collisions the two 
pion sources separated by a long distance are the direct and long-lived resonance 
sources. 

It is interesting to compare at the time evolutions of various observables in 
the context of the HRM for RHIC Au-|-Au collisions. Figure 18 shows the time 
evolution of the pion elliptic flow, pion HBT, and slope parameters from the 
rescattering calculation with an impact parameter of 8 fm and averaged over 100 
events. Lines fitted to the points from the calculation are shown for convenience. 
All quantities are extracted at midrapidity, elliptic flow includes all p-p, HBT is 
calculated for 125 < pr < 225 MeV/c, and the slope parameters are calculated in 
the region niT < 1 GeV. In these calculations, rescattering drives all of the time 
evolution seen in the various quantities, i.e. if rescattering were turned off, there 
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Fig. 14. Preezeout (and decay) time distributions for pions, kaons, and nucleons for SPS Pb+Pb 
(a) from a full rescattering calculation, and (b) with rescattering turned off (Prom Ref.[21]). 



would be no change in their t = values for t > 0. It is seen that the elliptic flow 
develops the earliest, stabilizing at about 5 fm/c, the HBT parameters stabilizing 
next at about 10 fm/c, and the latest being the slope parameters which require a 
time somewhat longer than 25 fm/c to stabilize (note that rescattering calculations 
are carried out to a time of 100 fm/c). 

3.3.2. Subdivision test of HRM code 

Before presenting comparisons of the HRM code with experimental results, it is 

worthwhile to address two criticisms which have been made against using the HRM 
to make RHIC-energy predictions. They are 1) the initial state for the calculation is 
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Fig. 15. Transverse mass distribution at freeze out from HRM for SPS Pb+Pb collisions for pions, 
kaons, and nucleons at midrapidity with exponential fit (Prom Ref.[21]). 




Fig. 16. Time evolution of the slope parameters in HRM for SPS Pb+Pb collisions for a full calcu- 
lation (solid lines) and for a "pion gas" calculation (dashed line). Experimental slope parameters 
from NA44 (points) are plotted at t = 70 fm/c for comparison (Prom Ref. [21]). 



too hot and too dense to be considering hadrons, thus the results are meaningless, 
and 2) the calculational results may have reasonable agreement with data but that 
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Fig. 17. Projections of tlic tliree-dimensional two-pion correlation function for SPS Pb+ Pb onto 
the QTside: QTouti Q Long 3^6S from HRM. A fit to Eq. (15) is also shown (Prom Ref.[21]). 



is only accidental because the calculation is dominated by computational artifacts 
which strongly influence the results. A response to 1) is that the HRM should be 
viewed as a limiting case study of how far one can get with an extreme and simple 

model such that maybe we can learn something about the true initial state from this 
unexpected agreement with data, e.g. maybe hadron-like objects exist in the QGP, 
and/or the QGP has a short lifetime and then quickly hadronizes. The response to 
2) is that the HRM has been tested for Boltzmann-transport-equation-like behav- 
ior and the influence of superluminal artifacts using the subdivision niethod.^^'"^'^ 
Although in that test the HRM results were shown to not be significantly affected 
by using a subdivision of Z = 5, it was not studied whether I = 5 was sufficiently 
large to signiflcantly reduce the superluminal artifacts to make the test meaning- 
ful. Superluminal artifacts can be introduced into a scattering code at the point at 
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Fig. 18. Time evolution of V2, HBT, and tut slope pajrameters from HRM for RHIC Au+Au. 

which the scattering cross section is used to determine whether two particles are 
sufficiently close in space and time for a scattering to take place if the density of 
particles in the calculation is sufficiently high.^^ In a subdivision test of a scattering 
code, the particle density, p, is increased by a factor I, the subdivision, while at the 
same time the scattering cross section, a, is decreased by this same factor, i.e. 



(22) 
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Running the scattering code in this configuration with Z > 1 in principle reduces the 
supcrhiminal artifacts and tests whether the code is properly solving the Boltzmann 
transport equation. 

In the present test, subdivisions oil = 1, 5, and 8 are used in HRM calculations 
of RHIC-energy Au+Au collisions with 6 = 8 fm centrality. Plots of the transverse 

signal velocity distributions for these subdivision are shown in order to determine 
how effective these subdivisions are in reducing the super luminal effects. The pion 
HBT radius parameters are also shown for these subdivisions in order to deter- 
mine how sensitively they depend on them. The left plot in Figure 19 shows the 
transverse signal velocity distributions for all particles in the calculation. As shown, 
there are indeed superluminal effects present for Z = 1, but the higher subdivisions 
significantly reduce these effects. Prom this, it is seen that using I = 5 and I = 8 
should each provide a valid test of the effects of these artifacts on the results of 
the calculation. Results from calculating pion HBT parameters vs. px for 6 = 8 fm 
is shown in the right plot of Figure 19. As is seen, the higher subdivisions do not 
significantly effect the HBT results. Radial and elliptic flow results from the HRM 
calciilations can also be shown to not be effected significantly by using these higher 
subdivisions.'*" Thus one can conclude from this test that the previously published 
results and present results from the HRM are not affected by superluminal artifacts, 
and criticism 2) above is answered. One can now proceed with comparisons with 
data below. 
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Fig. 19. Subdivision test of the HRM code, showing transverse signal velocity distribution and 
pion HBT parameters vs. pt for 6 = 8 fm centrality and subdivisions 1, 5, and 8. 
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3.4. Comparison of HRM with AGS (y/s ~ 2 GeV/nucleon) and 
SPS (y/s ~ 20 GeV/nucleon) experiments 

Having discussed some of the features of the HRM, comparisons between HBT 
results from the HRM with AGS and SPS experiments are now presented below in 
Figures 20-23. 

Figure 20 shows a comparison between boson source parameters from the HRM 
(triangles) with those from the fixed-target central collisions 200 GcV/nucleon 
S-|-Pb measured by SPS experiment NA44^^''^^ and 14.6 GeV/nucleon Si+Au from 
the AGS experiment E859/E866^^ (circles). Both experiments are based on small 
acceptance magnetic spectrometers with good particle identification allowing them 
both to also carry out HBT measurements with kaon pairs as well as pion pairs as 
shown in the figure. As seen, NA44 also measured pion source parameters with a 
low-pT cut (< pt >'~ 150 MeV/c) and a high-pT cut (< pr >'~ 450 MeV/c) on the 
pion momentum. The HRM is seen to follow the trends of the data rather well, pre- 
dicting the decrease in the source parameters with the higher pion momentum cut 
(NA44) and for both experiments predicting smaller source parameters extracted 
from the kaon-pair HBT measurements. In the context of the HRM, the explanation 
for the smaller parameter sizes extracted used kaon pairs is due to the generally 
smaller scattering cross sections iov K — n and K — N reactions compared with 
those for pions. The dashed lines in the figure show the parameters extracted from 
HRM with rescattering turned off, showing that rescattering plays a crucial role in 
determining the scale and kinematic dependencies of the boson source parameters 
measured by HBT interferometry. 

Figure 21 presents a comparison between the HRM for various model configura- 
tions and HBT parameters extracted from fixed-target central 158 GeV/nucleon 
Pb+Pb collisions from SPS experiment NA49 (based on a large acceptance 
detector). '^'^ The dashed lines are projections of the NA49 data points to guide the 
eye. For the purposes of this figure, "IOC" refers to the use of Eq. (20) in the HRM 
(the normal method used) and "pill" refers to taking thad = and z^ad uniformly 
distributed in the region z = ±1. As seen, the HRM predictions show robustness to 
the first three running configurations predictions being close to the measurements 
for the cases "IOC", "pill", and "pill(old)", the latter being a calculation without 
resonances which accounts for the A parameter being unity for that case. The cases 
"IOC (pion gas)" and "IOC (no RS)" , referring to running the HRM code without 
kaons and nucleons and with rescattering turned off, respectively, are seen to de- 
viate significantly from the measurements showing the importance of including the 
pions and kaons and rescattering in the calculations. 

Figures 22 and 23 show more comparisons of HBT source parameters extracted 
from the HRM and the SPS NA49 and NA44 experiments. Figure 22 shows a com- 
parison between HRM (points) and the trends of NA49 central 158 GeV/nucleon 
Pb-I-Pb experimental results^'^ (dashed lines with same meaning as in Figure 21) 
for the kr dependence of the pion source parameters. Figure 23 compares an 
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Fig. 20. Comparison between boson source parameters from the HRM (triangles) with those from 
SPS experiment NA44 and AGS experiment E859/E866 (circles). The dashed lines show the 
parameters extracted from HRM with rescattering turned off (From Ref. [26]). 



overview of HBT-cxtractcd pion and kaon source parameters measured in the NA44 
experiment^^'^"*'^^'^^ with those calculated with the HRM. Although there are some 
minor disagreements, the HRM is seen overall to qualitatively describe the trends 
of the data for both experiments. 



3.5. Comparison of HRM with RHIC ~ 130 - 200 

GeV/nucleon) experiments 

Results from the first year of running of the Relativistic Heavy Ion Collider (RHIC) 
for Au+Au collisions at ^/s = 130 GeV showed surprisingly large pion elliptic flow'^^ 
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Fig. 21. Comparison between SPS experiment NA49 HBT results for Pb + Pb with HRM for 
various model configurations. The dashed lines are projections of the NA49 data points to guide 
the eye (From Ref.[21]). 
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Fig. 22. Comparison between HRM (points) and the trends of NA49 Pb-i- Pb experimental results 
(dashed lines with same meaning as in Figure 21) for the fcy dependence of the pion source 
parameters (From Ref. [21]). 



and surprisingly small radii from two-pion HBT interferometry. ^"'^^ Hydrodynam- 
ical models agreed with the large elliptic flow seen in the RHIC data"^^ but signifi- 
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Fig. 23. Comparison between HRM calculations and boson source parameters extracted from the 
NA44 experiment (Prom Ref.[36]). 



cantly disagreed with the experimental HBT radii. On the other hand, relativis- 
tic quantum molecular dynamics calculations which include hadronic rescattering, 
for example RQMD v2.4,''^ significantly under-predict the elliptic flow seen in the 
RHIC data^^ but predict pion HBT radii comparable to the data.^^ A calcula- 
tion was made to extract HBT radii with a hydrodynamical model coupled with 
a hadronic rescattering afterburner with the result that the HBT radii were sig- 
nificantly larger than measurements.*^ This situation lead to the first big mystery 
from RHIC, sometimes called the "HBT Puzzle." It has even been suggested that 
one should call into question the current understanding of what information pion 
HBT measurements give.'*^ In this context, the HRM model is compared with both 
elliptic flow and HBT measurements from RHIC to see if it can shed any light on 
this situation. 

3.5.1. Elliptic flow results 

The elliptic flow variable for a collision, V2, is defined as^'' 

V2 =< cos(2(/>) > (23) 
(j) = arctan(— ) 

Px 

where and Py are the x and y components of the particle momentum, and x 
is in the impact parameter direction, i.e. reaction plane direction, and y is in the 
direction perpendicular to the reaction plane. 
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Fig. 24. Calculations of V2 from the HRM for 6 = 8 fm for pions and nucleons compared with 
STAR measurements at 11-45% centrality. The plotted points with error bars are the HRM cal- 
culations and the lines show the trends of the STAR measurements. Average errors on the STAR 
measurements are < 0.002 for pions and 0.006 for protons-|-antiprotons (Prom Ref.[22]). 



Figure 24 shows the pT dependence of V2 for pions and nucleons extracted from 
the 5 = 8 fm HRM calculation compared with the trends of the STAR measurements 
for 7r+ -|- TT^ and p + p a.t 11 — 45% centrality, "^^ which roughly corresponds to this 
impact parameter. 

Figure 25 compares the pr dependence of V2 for kaons from the 5 = 8 fm 
HRM calculation with the STAR measurements for K^g at 11 — 45% centrality.'*^ 
As seen, the HRM calculation values are in reasonable agreement with the STAR 
measurements. The flattening out of the pion and nucleon V2 distributions for 
Pt > 2 GeV/c is consistent with that seen in STAR and PHENIX results for 
minimum-bias hadrons^'''^* (the kaon V2 calculation does not extend higher than 2 
GeV/c in Figure 5 due to limited statistics). Thus, the same rescattering mechanism 
that can account for the radial flow seen in rri'r distributions, e.g. Figures 15 and 
16, also is seen to account for the magnitude and pt dependence of the elliptic flow 
for pions, kaons, and nucleons. 



3.5.2. "Year- 1 " HBT results 

The pion source parameters extracted from HBT analyses of HRM calculations 

for three different impact parameters, b = 0, 5, and 8 fm, are compared with 
STAR 7r~ measurements at three centrality bins^° in Figure 26. Note that the 
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Fig. 25. Calculation of V2 from HRM for kaons at 6 = 8 fm compared with STAR measurements 
for at 11-45% centrality (Prom Ref.[22]). 



PHENIX experiment HBT results"^* are in basic agreement with the STAR results. 
The STAR centraUty bins labeled "3", "2", and "1" in the figure correspond to 12% 
of central, the next 20%, and the next 40%, respectively. These bins are roughly 
approximated by the impact parameters used in the HRM calculations, i.e. the 
average impact parameters of the STAR centrality bins are estimated to be within 
±2 fm of the HRM calculation impact parameters used to compare with them. In 
the left panel, the centrality dependence of the HBT parameters is plotted for a 
Pt bin of 0.125 — 0.225 GeV/c. In the right panel, the tot dependence of the HBT 
parameters is plotted for centrality bin 3, for the STAR measurements, or 6 = 
fm, for the HRM calculations. Although there are differences in some of the details, 
the trends of the STAR HBT measurements are seen to be described rather well 
by the HRM calculation. 



3.5.3. Azimuthal HBT results 

It was first shown experimentally at the AGS that looking at the azimuthal de- 
pendence of HBT with respect to the reaction plane can give another handle on 
the space-time evolution of the pion source.*^ For central c;ollisions with vanishing 
impact parameter (6 ~ 0) one would expect no azimuthal dependence since the 
initial collision is symmetric about the beam axis. However, for non-central {b > 0) 
collisions, a definite initial asymmetry with respect to the beam axis will exist which 
may be reflected in the extracted HBT source parameters. The STAR collaboration 
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Fig. 26. Comparison of HBT source pcirameters from HRM with STAR measurements as a func- 
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the rescattering calculations are the open symbols. The errors on the STAR measurements are 
statistical+systematic (Prom Ref.[22]). 



has recently published experimental results on the azimuthal dependence of pion 
HBT parameters in central and non-central RHIC collisions. ^'^ Preliminary calcu- 
lations have been made with the HRM to extract azimuthal HBT parameters to 
compare with the STAR results. These comparisons are shown in Figures 27 and 28 
below. In Figure 27 the azimuthal dependence of HBT parameters is shown for 6 = 
fm rescattering model calculations compared with results from STAR central colli- 
sions (0 — 5% centrality). As seen for both the calculation and data, no oscillations 
occur with respect to (j) for any of the parameters except the cross-term parameter, 
^outside^ as would be expected for an azimuthally symmetric system. Figure 28 
shows a similar comparison for 6 = 4 fm calculations and STAR medium-central 
collisions (10 — 20% centrality). As would be expected, the non-central collisions 
break the azimuthal symmetry of the pion source and oscillations are now seen in 
all parameters. The calculations are seen to be in reasonable qualitative agreement 
with the data. 



3.5.4. Discussion of RHIC results 

As shown above, the elliptic flow as well as the features of the HBT measurements at 
RHIC can be adequately described by HRM with the hadronization model parame- 
ters given earlier in Table I. The results of the calculations arc found to be sensitive 
to the value of Thaj. used, as was studied in detail for SPS HRM calculations.^^ For 
calculations with Thad > 1 fm/c the initial hadron density is smaller, fewer collisions 
occur, and the rescattering-generated flow is reduced, reducing in magnitude the 
radial and elliptic flow and most of the HBT obscrvables. Only the HBT parameter 
Rhong increases for larger Thad reflecting the increased longitudinal size of the initial 
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Fig. 27. Azimuthal pion HBT parameters from the rescattering model with 6 = fm centrahty 
compared with central STAR results 



hadron source, as seen in Eq. (20). One can compensate for this reduced flow in the 
other observables by introducing an ad hoc initial "flow velocity parameter" , but 
the increased Rhong cannot be compensated by this new parameter. In this sense, 
the initial hadron model used in the present calculations with T^a^ = 1 fm/c and 
no initial flow is uniquely determined with the help of Rhong- 

At this point, one can consider the physical significance of the HRM in two dif- 
ferent ways. The first way is to accept that it is physically valid in the time range 
where hadronic rescattering should be valid, e.g. for times later than when the parti- 
cle density reaches about 1 fm.~^, and to take the initial state hadronization model 
as merely a parameterization useful to fit the data. Considering the calculation this 
way, one can at least expect to gain an insight into the phase-space configuration 
of the system relatively early in the collision (in the calculation 1 fmr'^ occurs at a 
time 4 fm/c after hadronization). The second way to consider the calculation is to 
see if it is possible to also physically motivate the initial state hadronization model. 
An attempt to do this is given below. 

In order to consider the present initial state model as a physical picture, one 
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Fig. 28. Azimuthal pion HBT parameters from the rescattering model with 6 = 4 fm centrality 
compared with medium-central STAR results 



must assume: 1) hadronization occurs very rapidly after the nuclei have passed 
through each other, i.e. Thad = 1 fm/c, 2) hadrons or at least hadron-like objects 
can exist in the early stage of the collision where the maximum value of p approaches 
8 GeY/fm^, and 3) the initial kinetic energies of hadrons can be large enough to 
be described by T = 300 MeV in Equation 1. 

Addressing assumption 2) first, in the calculation the maximum number density 
of hadrons at mid-rapidity at t = fm/c is 6.8 fm~^, rapidly dropping to about 
1 fm~^ at t = 4 fm/c. Since most of these hadrons are pions, it is useful as a 
comparison to estimate the effective volume of a pion in the context of the n — n 
scattering cross section, which is about 0.8 fm^ for s-waves.^* The "radius" of a pion 
is found to be 0.25 fm and the effective pion volume is 0.065 fm^, the reciprocal 
of which is about 15 fmr'^. From this it is seen that at the maximum hadron 
number density in the calculation, the particle occupancy of space is estimated to 
be less than 50%, falling rapidly with time. One could speculate that this may be 
enough spacial separation to allow individual hadrons or hadron-like objects to keep 
their identities and not melt into quark matter, resulting in a "super-heated" semi- 
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classical gas of hadrons at very early times, as assumed in the present calculation. 

Since the calculation takes the point of view of being purely hadronic, it is in- 
structive to consider assumptions 1) and 3) in the context of the Hagedorn thermo- 
dynamic model of hadronic collisions. According to Hagedorn, the mass spectrum 
of hadrons of mass m increases proportional to exp(m/To) in hadronic collisions, 
where Tq = 160 MeV is the limiting temperature of the system. This seems to con- 
tradict the value T = 300 MeV needed in the present case in Eq. (18) to describe 
the data. The Hagedorn model assumes that a) the system comes to equilibrium 
and b) the details of particle production via direct processes and through resonance 
decay average out. Neither of these assumptions is necessarily guaranteed at very 
early times in the collision. The use of the thermal functional form, Equation 1, 
to set up the initial transverse momenta of the hadrons in the present calculations 
is convenient but not required. For example, the exponential form exp (— mT/T'e) 
(where Te is a slope parameter) which docs not describe thermal equilibrium, could 
have equally well been used. This exponential form of the transverse mass distri- 
bution was successfully used previously in rescattering calculations to describe SPS 
data.26>27 

Assumption 1) can also be motivated by the Color Glass Condensate model. ^^'^^ 
In the usual version of this picture, after the collision takes place the Color Glass 
melts into quarks and gluons in a timcscalc of about 0.3 fm/c at RHIC energy, and 
then the matter expands and thcrmalizcs into quark matter by about 1 fm/c. In 
the context of the HRM calculations, it is tempting to modify the collision scenario 
such that instead of the Color Glass melting into quarks and gluons just after the 
collision, the sudden impact of the collision "shatters" it directly into hadronic 
fragments on the same timescale as in the parton scenario due to the hadronic 
strong interactions. 

3.6. Predictions from HRM for LHC (y/s ~ 5500 GeV/nucleon) 
experiments 

Since it has been found that the predictions from HRM agree rather well with AGS, 
SPS, and RHIC measurements, it is interesting to use this model to make similar 
predictions for Pb+Pb collisions at the LHC. Preliminary calculations for the LHC 
have been carried out with HRM, the results from which are shown below. In per- 
forming LHC calculations, the following parameters were used in the code: 1) a 
collision impact parameter of 6 = 8 fm, 2) an initial temperature parameter of 500 
MeV, 3) a hadronization proper time for the initial system of 1 fm/c, 4) a dn/dy at 
mid-rapidity for central collisions for all particles of 4000, and 5) an initial rapidity 
width of 4.2. These parameters were judged to be reasonable guesses to simiilate 
LHC Pb-I-Pb collisions. They at least satisfy the self-consistency check that sum- 
ming over the energy of all particles in an event at the end of the calculation agrees 
with the input total energy of a LHC Pb+Pb collision with an impact parameter of 
6 = 8 fm. An impact parameter of 6 = 8 fm was chosen for the present preliminary 
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study both to obtain non-negligible elliptic flow values and for calculational conve- 
nience (even for this impact parameter the cpu time used by the code for each LHC 
Pb+Pb event was about 60 hours). For item 3) above, the hadronization proper 
time was taken to be the same as was used in the SPS and RHIC calculations. 
Results of these calculations are compared with similar calculations at 6 = 8 fm 
centrality for RHIC Au-|-Au collisions and are shown in Figures 29, 30, and 31. All 
of these results are obtained at mid-rapidity, i.e — 2 < y < 2. 

Figure 29 shows the pion freezeout time and z-position distributions for LHC 
Pb-|-Pb and RHIC Au-I-Au from HRM. As seen, the average pion freezeout time 
and z-position for LHC are about twice as large as those for RHIC. Although the 
tails of the freezeout time distributions extend beyond 100 fm/c, the peaks for the 
LHC and RHIC occur at fairly short times in the collision, at about 5 fm/c and 
10 fm/c, respectively. Thus, efi^ects from earlier times in general have the greatest 
influence on the results from these calculations. 

Figure 30 shows the radial and elliptic flow predictions for LHC Pb+Pb com- 
pared with RHIC Au-|-Au from HRM. As seen in the mr-distribution plot on the 
left, although all species of particles start from a common temperature in the cal- 
culation, after rescattcring the exponential slope parameters follow the usual radial 
flow pattern of slope(v:) < slope{K) < slope{N) for both LHC and RHIC. The 
slopes are seen to be consistently larger at LHC than RHIC, as well as for the 
degree of radial flow which is built up. Looking at the plot of pion elliptic flow vs. 
Pt on the right, it is seen that LHC and RHIC give about the same values. This is 
due to the elliptic flow stabilizing at a very early stage in the HRM calculation. 

Figure 31 shows the pion HBT parameters vs. pT for LHC Pb+Pb compared 
with RHIC Au+Au from HRM. The transverse radius parameters, Rrside and 
Rxout, are seen to be somewhat larger and show a stronger px dependence for 
LHC as compared with RHIC. The longitudinal radius parameter, Rhong is seen to 
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Fig. 29. Pion freezout time and z-position distributions for LHC Pb+Pb and RHIC Au+Au for 
6 = 8 fm centrality collisions at mid-rapidity from HRM. 
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Fig. 30. rriT distribution and pion elliptic flow predictions for LHC Pb+Pb compajred with RHIC 
Au+Au for 6 = 8 fm centrality collisions at mid-rapidity from HRM. Note that in the plot on the 
left, corresponding RHIC slope parameters from HRM are show in parentheses for comparison 
with the LHC values. 

be significantly larger for LHC as compared with RHIC, clearly reflecting that the 
pion freezeout times at LHC are twice as long as at RHIC according to HRM. The 

A parameter is seen to increase with increasing pT in the same way for both LHC 
and RHIC, reflected the reduced influence of long-lived resonances at the higher pt 
values. 

Summarizing the results of this preliminary study, it is predicted from HRM 
that medium-peripheral (6 = 8 fm) LHC Pb-|-Pb collisions will produce more radial 
flow and larger HBT radii than the analogous RHIC Au-|-Au collisions, although 
elliptic flow and the A parameter values will look the same. 

4. Summary 

Over the past 25 years the HBT interferometry technique as applied to relativistic 
heavy ion collisions has evolved from being little more than a curiosity to becoming 
a standard tool applied by most RHI experiments to help understand their data. It 
has evolved in both its level of theoretical sophistication as well as in the level of 
precision at which HBT measurements can be carried out. The issue still remains 
of how to physically interpret the HBT observables which are measured. As is the 
clear theme of the present review, hadronic scattering models provide a means of 
addressing this issue of interpretation. Because of their success in describing the be- 
havior of experimental HBT parameters, in addition to other measured observables 
such as the radial and elliptic flow, such models can be used in conjunction with 
experiments to view relativistic heavy ion collision at an earlier stage prior to the 
randomizing effects of the rescattering process. It will be more than interesting to 
carry out HBT measurements at the LHC in several years and thus to see if nature 
presents us with a host of new HBT puzzles to solve. 
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